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Abstract 

We consider the standard model extended by heavy right handed fermions transforming as 
triplets under S\J{2)l, which generate neutrino masses through the Type-III seesaw mechanism. 
At energies below their respective mass scales, the heavy fields get sequentially decoupled to give 
an effective dimension-5 operator. Above their mass thresholds, these fields also participate in 
the renormalization of the wavefunctions, masses and coupling constants. We compute the renor- 
malization group evolution of the effective neutrino mass matrix in this model, with particular 
emphasis on the threshold effects. The evolution equations are obtained in a basis of neutrino pa- 
rameters where all the quantities are well-defined everywhere, including at ^13 = 0. We also point 
out the important role of the threshold effects and Majorana phases in the evolution of mixing 
angles through illustrative examples. 
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I. INTRODUCTION 



In the last decade, results from solar, atmospheric, accelerator and reactor experiments 
looking for neutrino flavour oscillations have succeeded in establishing that atleast two of 
the neutrinos are massive and there is mixing between different flavors 1|. The present 
best-fit values of the mass squared differences and mixing angles determined from analyses 
of global data on neutrino oscillation are 



Ami, = 7-Q5-om X 10^^ eV^, |Amii| = 2.401J]:! x lO^^ eV^ 
in^ 9u = O.SOtrol , sin' ^23 = 0.501°:^! , sin^ ^^^3 = O.Oll°:°f 



where Amfj = mf — are the mass squared differences and 6ij the mixing angles. The 
relative position of the third mass eigenstate ms with respect to the other two is unknown, 
though the solar neutrino data give Arri^i > 0. This results in two possible orderings of the 
neutrino masses: normal (mi < m2 < m^) and inverted {m^ < nii < 1112). 

One of the most distinctive features emerging out of the above results is the occurrence of 
two large and one small mixing angles which is rather different from the quark sector where all 
three mixing angles are small. The absolute masses of neutrinos are also orders of magnitude 
smaller than those of quarks and charged leptons, the current bound from cosmology on the 
sum of neutrino masses being ^ mj ^ 1.5 eV [sl. The most favored mechanisms to generate 
such small neutrino masses and nontrivial mixings are the so called seesaw mechanisms which 
need the introduction of one or more heavy fields. At energies below their mass scales, the 
heavy fields get integrated out giving rise to an effective dimension-5 operator ^ 

= k^IlIl^xJ) , (1) 

where II and cj) are respectively the lepton and Higgs doublets belonging to the standard 
model (SM). Here ^5 is the effective coupling which has inverse mass dimension and can be 
expressed in terms of a dimensionless coupling as K5 = CI5/A, with A some high energy 
scale. In this picture the SM serves as an effective theory valid upto the mass scale A, 
which can be taken to be the mass of the lightest of the heavy fields. Such an operator 
violates lepton number by two units and hence gives rise to Majorana masses for neutrinos: 
mi, ~ Ks'y^, where v is the vacuum expectation value of the Higgs field after spontaneous 
symmetry breaking. Taking v ~ 246 GeV, a neutrino mass of ~ 0.05 eV implies A ~ 10^^ 
GeV if as ~ 1. 



There are four possible ways to form a dimension-5 gauge singlet term out of the two 
lepton doublets and two Higgs doublets: (i) each Il-4> pah forms a fermion singlet, (ii) each 
of the Il-^l and 0-0 pair forms a scalar triplet, (iii) each /x-0 pair forms a fermion triplet, 
and (iv) each of the Il-^l and 0-0 pair forms a scalar singlet. Case (i) can arise from the 
tree level exchange of a right handed fermion singlet and this corresponds to the Type-I 
seesaw mechanism Case (ii) arises when the heavy particle is a Higgs triplet giving 
rise to the Type-II seesaw mechanism 

aa. 

For case (iii) the exchanged particle should 
be a right-handed fermion triplet, which corresponds to generating neutrino mass through 
the Type-III seesaw mechanism 8|. The last scenario gives terms of the form i'iCl which 
cannot generate a neutrino mass. 

Type-III seesaw mechanism mediated by heavy fermion triplets transforming in the ad- 
joint representation has been considered earlier in j^, 9]. Very recently there has been a 
renewed interest in these type of models. The smallness of neutrino masses usually implies 
the mass of the heavy particle to be high ~ 10^^^^^ GeV. However, it is also possible to 
assume that one or more of the triplets have inasses near the TeV scale, making it possible 
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In such models, the Yukawa 



to search for their signatures at the LHC lol . 
couplings need to be small to suppress the neutrino mass. Lepton flavour violating decays 
in the context of Type-III seesaw models have also been considered in [l^. Recently it has 
also been suggested that the neutral member of the triplet can serve as the dark matter and 
can be instrumental in generating small neutrino mass radiatively jisl . 

The possibility of being able to add one triplet fermion per family without creating 
anomalies was one of the consequences of a general analysis in IGj which discussed adding 
an extra U(l) gauge group to the SM. Possible ways of adding fermion triplets in an anomaly- 



free manner have been explored [17|, with some specific models studied in [18|. A possible 
origin of such an extra U(l) gauge group has been proposed in [l^. Fermions in the adjoint 
representation fit naturally into the 24-dimensional representation of SU(5), and can rectify 
the two main problems encountered in SU(5) Grand Unified Theory (GUT) models, viz. 
generation of neutrino masses and gauge coupling unification 1^, [ill, UQ, 21]. The latter 
requirement constrains the fermionic triplets to be of mass below TeV for Mqut ~ 10^^ 



GeV, making the model testable at the LHC. Leptogenesis mediated by triplet fermions has 



been explored in 



22l |. Additional fermions transforming as triplet representations in the 



context of left-right symmetric model have been studied in 23(]. Minimal supersymmetric 
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standard model extended by triplet fermions has recently been considered in 24 1. 

Whether the exchanged particle at the high scale is a singlet fermion (Type-I seesaw) or 
a triplet fermion (Type-Ill seesaw), the light neutrino mass matrix is given as m^M^^mo. 
Here m/) is the Dirac mass matrix coupling the left handed neutrinos with the right handed 
heavy fields, and M.^ is the Majorana mass matrix for the right handed fields. Thus the 
generation of the light neutrino mass matrix is similar in the Type-I and Type-Ill seesaw 
mechanisms, both of which are fermion mediated. Since the neutrino mass is generated at 
the high scale while the neutrino masses and mixings are measured experimentally at a low 
scale, the renormalization group (RG) evolution effects need to be included. These radiative 
corrections in Type-I and Type-Ill seesaw are different, since the heavy fermions couple 
differently to the other particles in the theory. We note that below the mass scale of the 
lightest of the heavy particles, the effect of all heavy degrees of freedom are integrated out 
and the effective mass operators in these scenarios become identical. 

The effect of RG induced quantum corrections on leptonic masses and mixings have been 



studied extensively in the literature [25 
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32|. 



interesting consequences such as the generation o: 



small mass splittings for degenerate neutrinos 
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ar 



m 



e mixing angles 
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41 



generation of ^13 starting from a zero value at the high scale 
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hese effects can have 
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49|, 



34 



LJ46 



50, 



35, 



36, 



3d, 



38|, 



47] , or radiative 
5l|. RG induced 



deviations from various high scale symmetries like tri-bimaximal mixing scenario 



or quark-lepton complimentarity 54, 



55 



56, 



52,1531,15 



571] and correlations with low scale observables 



have been explored. Such effects can have significant contributions from the threshold 
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60l |. The RG evolution of the neutrino mass operator in the SM and 



corrections 

;he Minimal Supersymmetric Standard Model (MSSM) in the context of Type-I seesaw 
59I, 
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62| and Type -II seesaw [63|, |6J] have been studied in the literature. In the context 
of Type-Ill seesaw with degenerate heavy fermions, the impact of the RG evolution on the 
vacuum stability and perturbativity bounds of the Higgs Boson has been explored in 65 1. 

In this work we study the RG evolution in the SM in the context of the Type-Ill seesaw 
model with nondegenerate heavy fermions. Our model consists of the SM with additional 
massive fermion triplets S with masses ~ Mj, [i = 1,2, ■■ ■ ,Mr) such that Mi < M2 < 
■ ■ ■ < Mr- Below the mass scale Mi all the triplets will be decoupled from the model and 
the RG evolution will be according to the SM. The triplets will manifest themselves at 



this low scale in the form of an effective operator k obtained by integrating out the heavy 
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fields. For energy scales above Mi, the effect of the heavy fermions will come into play 
successively and above Mj. all the three triplets will contribute to the RG running. We 
evaluate the contributions of these fermion triplets to the wavefunction, mass and coupling 
constant renormalization of the SM fields and of the triplet fields themselves. We obtain the 
/3-functions for RG evolution of the Yukawa couplings, the Higgs self-coupling, the Majorana 
mass matrix of the fermion triplets, the effective vertex n and the gauge couplings, including 
the extra contribution due to the additional triplets wherever applicable. We obtain analytic 
expressions for the runnings of the masses, mixing angles and phases in a basis where all the 



quantities are well-defined at every point in the parameter space including ^^13 = [51[. We 
also solve the RG equations numerically and present some illustrative examples of running 
of masses and mixing angles. We analyze the effect of the seesaw thresholds and Majorana 
phases and check if such a scheme can generate masses and mixing angles consistent with 
the current bounds. 

The plan of the paper is as follows. In Sec. [TTl we outline the basic features of the Type- 
Ill seesaw model including extra SU(2)2,-triplet fermions, and describe how the effective 
neutrino mass operator can arise by this mechanism. In Sec. IIIIl we describe how to include 
the varying mass thresholds of the heavy particles in the analysis, discuss the renormalization 
of the SM extended with heavy triplets, and give the expressions of the /3 functions including 
the effect of the extra triplets. In Sec. IIVI we detail the changes in the RG equations of the 
effective neutrino mass operator due to the inclusion of the extra fermion triplets. In Sec. |Vl 
we numerically demonstrate the modifications in the RG equations of the neutrino masses 
and mixing angles. We summarize our results in Sec. IVIi 

II. THE TYPE-III SEESAW MODEL 

We consider the Type-Ill seesaw model where three heavy fermions are added to each 
family of the SM. These fermions have zero weak hypercharge, i.e. they are singlets of the 
gauge group U(l)y of the SM. However, under the SU(2)i gauge, they transform as a triplet 
in the adjoint representation. In the basis of the Pauli matrices {a^, cr^, a^}, this triplet 
can be represented as 



-'R 



(2) 
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where = (E|j ^ iE|j)-\/2. For the sake of simphcity of further calculations, we combine 
Eij with its charge conjugate 

^ / j:f/V2 E^,^ \ ^ EgV 
and use the quantity E, defined as 

E = Efl + (Efl)^. (4) 

Clearly, E also transforms in the adjoint representation of SU(2) l. Note that though formally 
E = E^, the individual elements of E are not all Majorana particles. While the diagonal 
elements of E are indeed Majorana spinors which represent the neutral component of E, the 
off-diagonal elements are charged Dirac spinors. 



A. The Lagrangian 

Introduction of the fermionic triplets E will introduce new terms in the Lagrangian. The 
net Lagrangian is 

>C = CsM + >Cs , (5) 



where 



Here, 



where 



C-E — Cs^kin + ^'EjTnass + ^'E,Yukawa ■ (6) 



>CE,fem = Tr[S2^E] , (7) 
>Cs,ma.s = ^IV[EMsE] , (8) 



(10) 

is the completely anti-symmetric tensor in the S\J{2)l space. Here we have not written the 
generation indices explicitly. Ms is the Majorana mass matrix of the heavy fermion triplets 
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and Is is the Yukawa coupling. The SM fields II, and are S\J{2)l doublets and can be 
written as 

Each member of the SU{2)l doublet II is a 4-component Dirac spinor. Since the fermion 
triplet S is in the adjoint representation of SU(2)l, the covariant derivative of S is defined 

as 

D,^ = ^,J: + ^g2[W„^ , (12) 

where g2 is the SU(2) i gauge coupling. 

All the Feynman diagrams for the new vertices involving the triplet fermionic field E are 
given in the Appendix |X1 The Feynman diagrams for the SM particles are shown in the 
Appendix [Bl 



B. The effective vertex 



In the low energy limit of the extended standard model, we have an effective theory 
which will be described by the SM Lagrangian with the additional operators obtained by 
integrating out the heavy fermion triplets added to it. The lowest dimensional one of such 
operators is the dimension-5 operator^ 

= (W^cx'e<j)^ (0Vdi) + h.c., (13) 
= -'^fg {^Li4>Jib4>dj ^ {£ac£bd + £ab£cd) + h.c. , (14) 

where k is a symmetric complex matrix with mass dimension (—1). Generation indices 
f,gE {1, 2, 3} are shown explicitly and a, b,c,d E {1, 2} are the SU(2) ^ indices. In writing 
the last line we have used 

(0-*)af)(0"')cd = 25ad5bc " ^ab^cd 
=^ {(y'£)ba{(y'e)dc = ^Eda^bc " ^ba^dc (15) 

^ We use this form to emphasize the triplet nature of the pairs. Since all the dimension-5 operators 
are equivalent, we choose the normalization such that Kfg defined here matches that in ,66.] . 
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(a) 



FIG. 1: The effective vertex k at an energy <^ Mi, after all the heavy fermions have been 
decoupled from the theory. /, g G {1,2,3} are the generation indices. The SU(2)i and generation 
indices for S are not shown explicitly since they are summed over. 

and utilizing the (pd 4>a symmetry, we can write 

'^^da^bc — ^ba^dc = " {^ab^dc + ^db^acj ■ (16) 

The relevant diagrams in the complete theory giving rise to the effective operators in the 
low energy limit are shown in Fig[TJ The "shaded box" on the left hand side represents the 
effective low energy vertex k, while A(a) and A(^b) are the amplitudes of the diagrams labeled 
as (a) and (6) on the right hand side. The amplitudes are given by 

A(a) = ili' (If Ms^Fs)^^ [{e^cy')ab{e^cj%d] Pl , (17) 
A^b) = tfi' (FjM^^Fs)^^ [(^V)rf,(£V),,] Pl , (18) 

with e = A — D where D is the dimensionality that we introduce in order to use dimensional 
regularization. Note that A(^b) is obtained from A^a) just hj d ^ a interchange. Using 
Eq. f|T5|) one finally gets 

Aia) + A^b) = -il^' (yjM^^ls)^^ {eabScd + Sacebd) Pl ■ (19) 

This is equal to the left hand side of Fig. [1] with the identification 

K = 2yjMs^ls • (20) 

Equation (l20l) gives the Feynman rule for the low energy effective vertex k, as shown in 
the Appendix |Xl From Eqs. (I20l) and (fT^ . one gets the neutrino mass after spontaneous 
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symmetry breaking to be 



= - M^'Fs (21) 



which is the Type-III seesaw relation. Here, v denotes the vacuum expectation value of the 
Higgs field. 

III. RADIATIVE CORRECTIONS IN TYPE-III SEESAW 
A. Sequential decoupling of heavy fermions 

Let us consider the most general case when there are r triplets having masses Mi < M2 < 
■ ■ ■ < Mr-i < Mr- Above the heaviest mass M^, all the r-triplets are coupled to the theory 
and will contribute to the neutrino mass through seesaw mechanism as 

I .u 7)2 ('' + 1) (r + 1) (r + 1) 

mj = -- M^^ (/i > Mr) . (22) 

(r-+l) 

Here, is a [r x np] dimensional matrix (np is the number of fiavors, which is 3 in our 
case). Me is a [r x r] matrix and 101,^ is a [np x np] dimensional matrix. Below the scale 
Mr, the heaviest triplet decouples from the theory. Integrating out this degree of freedom 

(r) 

gives rise to an effective operator k. The matching condition at jj, = Mr is 



(r + l) 



(r + l) 



= 2Y^ (Mr)-' 

Mr 



(23) 

Mr 



This condition ensures the continuity of m,^ at = Mr. In order to get the value of the 
threshold Mr, we need to write the above matching condition in the basis where Ms = 
diag(Mi, M2, ■ ■ ■ , Mr). Here it is worth mentioning that the matching scale has to be found 
carefully since Ms itself runs with the energy scale, i.e. Mi = Mi{fi). The threshold scale 
Mi is therefore to be understood as Mj(yU = Mi). 

In the energy range Mr-i < /x < Mr, the effective mass of the neutrinos will be given as 

ml = - J ri^+2Yi M^' Y^ j . (24) 

The first term in Eq. fl2^ is the contribution of the integrated out triplet of mass Mr through 

(r) 

the effective operator k. The second term represents the contribution of the remaining (r — l) 
heavy fermion triplets, which are still coupled to the theory, through the seesaw mechanism. 

(r) (r) 

Ms is now a [(r — 1) x (r — 1)] matrix while Y-^ is a [(r — 1) x np] dimensional matrix. 



The matching condition at /i = Mr-i is 



(r-l) 

K 



(r) 

K 



Mr 



Mr 



(r) (r) 



(25) 



Generahzing the above sequence, we can say that if we consider the intermediate energy 
region between the {n — 1)**^ and the n*^ threshold, i.e. M„ > fi > M„_i, then all the heavy 
triplets from masses Mr down to M„ have been decoupled. In this region the Yukawa matrix 

(n) 

will be [(n — 1) x np] dimensional and will be given as 



ln—l,np 







V 



(n) 



(26) 



^ heavy triplets with masses 
M„ — Mr integrated out . 



(») 



Ms will be [{n — 1) x (n — 1)] dimensional. In this energy range the effective neutrino mass 
matrix will be 



(n) 

m. 



with 



(i) (")_, (n) (") 



(27) 



(28) 



while the matching condition at /i = M^ is given by Eq. (1251) with r replaced by (n + 1). 
For /i < Ml, all the heavy triplets will get decoupled and thus only W will contribute, which 
is the low energy effective neutrino mass operator. 



B. Dimensional regularization and renormalization 



Now we consider the radiative corrections to the fields, masses and couplings in our model, 
on the lines of that performed in 



29 



66l | in the context of Type-I seesaw. The wavefunction 



renormalizations are defined as 



fa 



(29) 
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where ifj G {IlAl^ ^r^ur, cIr} . We denote the renormahzed quantities as X and the corre- 
sponding bare fields as Xb- For the fermion triplets 



fg 



For the doublet Higgs 



whereas 



A 



B 



ZlA 



for the gauge bosons where A G {-B, W\ G"^}. For the Faddeev-Popov ghosts one has 

1 

cb = Z^c , 



(30) 

(31) 
(32) 

(33) 



however the ghosts will not appear in the RG evolution of the relevant quantities at one 
loop level. We introduce the abbreviation 



5Z 



X 



(34) 



where Zx denotes the renormalization constant of any of the relevant quantities X. 

We will use the dimensional regularization and the minimal subtraction scheme for renor- 
malization. In this renormalization formalism, the counter terms are defined such that they 
only cancel out the divergent parts. Thus the renormalization constants are of the form 

+ (35) 



k>l 



where the 5Zx,k are independent of e. In our scenario, at the one loop level, the renormal- 
ization constants are proportional to 1/e. The final results of course will be independent of 
the particular regularization as well as the renormalization scheme used for the calculations. 

The diagrams contributing to the renormalization constants of the different quantities are 
all shown explicitly in Appendix O The renormalization constants of different quantities 
are given by 



5Z^ 



167 



5Z, 



5Zs 



1 



T / ("),(") \ 1 



167r 



(36) 
(37) 
(38) 

(39) 
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where we have used the gauge, and the GUT normahzation of the gauge couphngs [25|. 
The Yukawa couphngs are renormahzed as^ 



1 / (") Q \ 1 



IGtt 
16712 [^^^ 



(40) 
(41) 



while the Majorana neutrino mass matrix gets renormahzed as 

5Zm, = --^(12 + 46)^72'-- 



(42) 



The addition of the right handed fermion triplets to the SM will contribute one extra diagram 
to the renormalization of the Higgs self-coupling A, as shown in the diagram (Gl) of the 
Appendix O This contribution will be^ 



5i 



47^2 



(n) 



Tr 



(n) (n) (n) (n) 

Yy%Y^% 



(43) 



Finally for the effective vertex k, the renormalization constant is 



OK - 



167r2 



2'"^{YX)+2{YXy'"^'-^^^ 



T(n) 



2 



- 6 Ul t - 



2 (") 



36 92 1^ 



1 

J e 



(44) 



We observe that there is no contribution from the fermion triplet S in the loop, which means 
that will not directly depend on the fermion triplets still coupled to the theory. However, 
during RG evolution an indirect dependence will creep in via the other couplings. 



C. Calculation of the (3 functions 

To calculate the f3 functions for the RG evolution of the Yukawa couplings, Majorana 
mass matrix, the effective vertex k and other relevant quantities, we consider the relations 

^ In jssi the contributions of fermion triplets to some of the above renormahzation constants are calculated 
in the context of SM extended with these fields. Their conventions of field normalizations are different 
and hence the results may differ upto numerical constants in certain cases. However, their Eq. (19) for 
SYi,, which is the same quantity as our SZy^ in Eq. (|¥T|l . is missing the Y^Ye term. The source of this 
term is the diagram labelled as {F2) in Ap pendix [Cl The extra contribution to SZye from the fermion 
triplets has also not been calculated in |65l| . 

Note that Ref. [gBi gives this quantity (^A in their Eq. (20)) to be of the form Tr(yjtFs). However, the 
additional contribution to the Higgs quartic coupling SZ\ should be of the form Ti{Y-^Yy,Y^Y^), since it 
comes from the diagram (Gl) in Appendix [Cl 
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between the bare (Xb) and the corresponding renormahzed (X) quantities given by 



yT 2 
1 

'2 



Z^^Yt^bZ^ Zi 
Zi^^ Z^kbZ^^Zi 



\ 1 

2 

L 

1 1_ 

2 
L 

1 1 
2 
L 



+ 5 k) 



(45) 
(46) 
(47) 
(48) 



where Z^;^ = PrZy.. We further use the functional differentiation method as in ^] to find 
the f3 functions for the Yukawa couphngs as 



1 c; (n) , (n) q q 



' ' \^ Tv^ I „^ „^ 



( -YX + -^vVs + T - -9? - - 



S ^2 



20 



9i 



92 



IQ-'Py. = Y^ [IyXu - IyX + T- il,? - - 8c 
16vr^/5y, = Y, ^^Y^ - h^u + T- \gl - \gl - 8,|) . 



(49) 
(50) 
(51) 
(52) 



Here 



T = Tr 



(n) , (n) 



YX + SrjFs + ^YXu + 3r]F, 



(53) 



(n) 



and j3x = iJi{dX / djj) . Note that Is is given in Eq. (126|1 . with (n — 1) the number of heavy 
fermion triplets still coupled to the theory. 

Since the fermion triplets have non-zero SU(2)i charge, they couple to the W bosons and 
hence will affect the RG evolution of the gauge coupling §2 via 



where 



19 A(n-1] 



h = -— + 



(54) 



(55) 



6 3 

Note that if the number of heavy fermion triplets is < 2, the value of 62 is always negative. 
On the other hand, if the number is > 3, then 62 becomes positive above the mass scale M3. 
Adding fermion triplets shifts the gi-g2 intersection to higher energy scales, and the g2-g3 
intersection to lower energy scales, as can be seen from Fig. [2l The exact situation would 
depend on the values of Mj. 
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FIG. 2: The solid (red) line and the dashed (green) lines show the energy scale variations of g\ and 
53 respectively in the SM, which is unaffected in Type-Ill seesaw. The dotted (blue) line gives the 
SM running of (72 j while dot-dashed (magenta), dot-dot-dashed (sky) and densely dotted (black) 
lines show the running if there were one, two or three fermion triplets respectively. 



The RG evolution of A is given by 

levr^/?, = - 3A 0(7? + 3(72') + 'igl + \ Q^?? + + 4AT 

-8 Tr[F;Fen^r, + 'iYlY^YlY^ + SFjr.rir,] - 20 T^Y^'SSk 



(56) 



As it is evident from Eq. fl56p . the last term is the new contribution to the /3-function from 
the heavy triplets still coupled to the theory. 

The RG evolution of the Major ana mass matrix of the heavy triplet fermions is given by 



167r^/?ME 



YyX^ Pl + Y^Y^ ) 



M, 



(n) (n) \ " / (n) (n) 

Y^Y^ Pl + IeFJ 1 Pn 



12(72Ms 



(57) 



where it is always possible to separate the components of different chirality to get the left- 
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chiral part as 



Y^Yy^ ) Ms + Ms ( Y^Yy^ ] - UgjUj: 



(58) 



since Pl + Pr = I. Thus all the /^-functions are gauge- independent, as they should be. The 
anomalous dimension of Ms is 



levr^ = My' 



(n) (n) 



Y^Y^ Pl + Y^Y^ ) Pn 



(n) (n) 



M, 



+ 



y^y^ ] pl+{ y^y^ i pr 



(59) 



Similar to the left-chiral component of (3^^ in Eq. (|58|) . the left-chiral component of is 
- IGtt^ Vm, = f FsFJJ Ms + ( ^^^e ) " 12(72' • (60) 



As seen from Eq. 0271) . the RG evolution of the light neutrino mass matrix in^ is controlled 



(n) 



by the evolutions of both k and Q, which are given by 



167r /3k = + P^ K + kP^ 



(71) (n) (n) 



with 



3 (") 3 , 

= ^yJy^ - ^yX 



a. 



2T + A - 3gl 



Q (n), (n) c; q q 

pq = ly.X + \yx ; «Q = 2T - - . 



(61) 
(62) 



(63) 



(64) 



IV. RG RUNNING OF NEUTRINO MASSES AND MIXING ANGLES 

To derive the RG evolution for the neutrino masses and mixings we follow the standard 
procedure 29|, |6l|- At any energy scale /i, the neutrino mass matrix -m.^ can be diagonalized 
by a unitary transformation via 



Uu{fi)^-m^{fi)U^{fi) = diag(mi(/i),m2(/i),m3(/i)) . 



In a basis where Yg is diagonal, the neutrino mixing matrix is given as 



(65) 



PMNS 



(66) 
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where ?7pmns is the Pontecorvo-Maki-Nakagawa-Sakata neutrino mixing matrix [67|, |68 |. 
From Eqs. fH9|) it is seen that above and between the thresholds, off-diagonal terms will be 
generated in even if we start with a diagonal at the high scale, due to the Y-^Y-^ terms. 
These terms will give additional contributions to the evolution of different parameters. In 
the presence of Y^ with off-diagonal entries, the neutrino mixing matrix will be given as 

t/pMNs = UlU, , (67) 

where Ue is the unitary matrix that diagonalizes Y^Y^, by a unitary transformation. t/pMNS 
is parameterized as [3, Q| 

?7pMNS = diag(e*^% e''% e''^) . U . ^mg{e-'^\ e-'^\ 1) , (68) 



with 



U 



^ C12C13 S12C13 Si3e ^ 

-C23S12 - S23Sl3Cl2e*'^ C23C12 - S23Sl3Si2e*'' S23C13 
y S23S12 - C23Si3Ci2e*'^ -S23C12 - C23Si3Si2e*'^ C23C13 J 



(69) 



Here Cij and Sij are the cosines and sines respectively of the mixing angle 9ij, 5 is the Dirac 
CP violating phase, (pi are the Majorana phases. The "flavor" phases 5e, 5^ and 5r do not 
play any role in the phenomenology of neutrino mixing. 

In this work, we consider r = 3 heavy fermion triplets, one for each generation. Then 1^ 
is a 3 X 3 matrix at high scale and is identically zero for /i < Mi. The RG evolution of the 
neutrino parameters is then controlled by 

167r2/3y, = Y,F + a,Y, , (70) 

167r^/?m, = P^'^u + + Oi^WCVy , (71) 

where 

P = CeY^Y, + Cj^YiXj, , (72) 

F = D,YX + Dj^Y^Yj, . (73) 

Eqs. (1701) and (17T]) are essentially the same as the /3-functions given in Eqs. (HO]) . and 
(!62|) . which we rewrite in the above form for later discussions. For /i > M3 and /i < Mi, the 
evolutions of Y^. and nij^ can be written in simple analytic forms, using Table [B Note that 
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Ce 
















/i > Ms 


5 
2 


3 
2 


3 
2 


15 
2 


r- 


- 


^fl2 
4^2 


9T - -^fl2 _ 9 2 
10^1 2^2 


/i < Ml 


3 
2 





3 
2 





T- 


4^1 


4^2 


2r + A - 3^1 



TABLE I: Coefficients of tlie /3-functions governing tlie running of neutrino masses and mixings in 
the energy regimes /i > M3 and /i < Mi. Tlie quantity T is defined in Eq. (I53p . 

for yU > M3 the running of the neutrino masses will be governed by I3q and so P in Eq. (1721) 
is the same as Pq as defined in Eq. ( l64l) . On the other hand, for /i < Mi, we have P = 
as given in Eq. fl63l) . P and P are 3x3 matrices, with the rows and columns representing 
generations. We denote the elements of P and F by Pfg and Ffg. The coefficient of Pfg and 
Ffg in the running of Y^. and can be read off directly from [611], since the structure of 
Eqs. ( ffOl) and (TfTl) remain the same both in Type-I and Type-III seesaw. The values of Pfg 
and Ffg themselves will however be different because of different underlying theories. The 
values of the relevant coefficients in Type-III seesaw are shown in Table [H 

If we consider the running equations in the basis V5 = {nii; 6*12, 6*13, 6*23; 4>i',S}, then both 



6 and 6 become ill-defined at 6*13 = 



29 



30| and as a consequence, 6'i3 also becomes ill- 



defined because of its 6 dependence. This is only an apparent singularity. One can get 



30| or by using the basis 



rid of it by imposing a particular value of cot 5 at 6'i3 = 
Vj = {nii] 9i2, 023, 9l^] (pi] Jcp, Jqf}: where the singularity does not appear at all [5l|. Here 
Jcp and J'qp are defined as 

-^CP = ^Si2Ci2S23C23Si3C^3Sin5 , (74) 
-^CP = ^■5l2Cl2S23C23Sl3Ci3COS(5 . (75) 

In the limit ^13 —>■ 0, Jcp, Jqp —>■ 0. From the point of view of the experiments also, the 
Jarlskog invariant Jcp is the quantity which appears in the probability expressions for CP 
violation in neutrino oscillation experiments, and is therefore directly measurable. J^p is 
needed in order to have complete information on S, since Jqp has no information on the sign 
of cos 5. We also choose to write the RG evolution for Of^ instead of 6*13 as is traditionally 
done. This quantity turns out to have a smooth behaviour at ^^13 = 0. Moreover, since 
^13 > by convention, the complete information about 6'i3 lies within Of^. The information 
about the Dirac phase will be present in Jqp, J'cp- 
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327r2 012 647r2 Bf^ 3271^ 



723 



Pii sin 2012 

P22 -Qt2 sin 20120^3 (^Ats - Ats^ sin 2^12 sin 2^23 (Q^3C?2 + Qt-A2} sin 2^23 

^^33 -QI2 sin 2^12 sig -(^^^3 -^^3) sin 2012 sin 2^23 -(Q^3C?2 + Q^3S?2) sin 2023 

ReP21 2Q+ cos 2012C23 4 (l+3cf2 + ^^3S?2) «23 (Q^g " Q^s) sin 2012^23 

RePsi -2Q+ COS2012S23 4(^+022 + ^+522)023 (Q^3 " Q^s) sin 20i2C23 
ReP32 Q^2 sin 2012 sin 2023 2 (1+ - ^+3) sin 20i2 cos 2023 2(3+022 + 2^35^2 ) cos 2023 
Im P21 45i2C23 4 (,Sf3cf 2 + -623512) S23 2 (^23 - ^is) sin 20i2S23 
ImPgi -45i2S23 4(5^3022+^23522)023 2(523- 5i3) sin 2012023 
ImP32 2(5^3-^^3) sin 2012 4 (5230^2 + ^13^12) 

TABLE II: Coefficients of Pjg in the RG evolution equations of the mixing angles 0i2, 0i3 and 023, 
in the limit 0i3 0. 



The expressions 



'or the running of masses and Majorana phases are the same as the 



ones obtained in 6II] for the Type-I seesaw mechanism. (See Tables 5, 6, and 14 therein. 



Note that 0j in our paper corresponds to in [6l|-) The running of masses and the 
Majorana phases does not depend on the Dirac phase to the lowest order in 6*13. Hence the 
RG evolution equations do not change with the change in basis Vs ^ Vj. Running of the 
two large mixing angles ^12 and ^23, as given in Table [Til is also the same as that in the Vs 
basis since the quantities Sij and Q^, defined as 

^ _ |m3±mie2*'?^i|2 ^ _ |m3 ± m2e2*'^2 12 ^ _ |m2e2*'?^2 mie2*'?^i |2 

= Aml,^ (1 + C) ' = Aml^ ' = A^, '^^^^ 

_ mim3sin2(/)i _m2m3sin202 _ mimg sin (20i - 2^2) 

'^--Am2,^(l + C)' - Aml^ ' A^, ' ^''^ 

depend on the mass eigenvalues and Majorana phases only. However the running of 6'23, as 
seen from the Table UTt depends on the quantities Afj, Sf^ defined as 

T± _ 4 {m\ + ml) J^p ± Smim^iJ'cp cos 20i + Jqp sin 20i) 



r± _ 4 {mj + m|) J^p ± 8m2m3( J^p cos 202 + ^cp sin 202) /-qv 
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Pn 








-45i2 cos 26*12 


P22 


-Aag+ 


20(^(7 - 2g-) 


4512C90 cos 26*12 


P33 


Aag- 


-20(^0" - 2Q-) 


451259^ COS 2^12 


ReP2i 


4S23^^ 


2s23(^o" + 2ac+) 


— 85i2C23 COS 2^12 cot 26*12 


RePsi 


4023^+ 


2c23(^r]" +2^+) 


85i2S23 COS 2^12 cot 2^12 


ReP32 


-2 sin 2^12 cos 2^23 


sin 2^12 cos 2^23 (^o" -2^^-) 


— 45i2 COS 2^12 sin 2^23 


ImPai 


2s23{g^ - 20+) 


4^23^+ 


-42^2^23 cot 26*12 


ImPai 


2023(^0+ -2g+) 


4c23a+ 


4Qr2'S23 cot 26112 


ImP32 


sin 26*12 (^o" +2^-) 


-2 sin 2^12^7 






TABLE III: Coefficients of Pfg in tlie RG evolution equations of tlie Jarlskog invariant JcP; t^ic 
quantity J^p = Jcp cot (5, and tlie Majorana phase difference — (/>2), in the limit ^13 0. The 
convention used here is a = S12C12S23C23, and Jcp = (a/2)si3Ci3 sin5. 

~± _ 4 {ml + ml) Jcp ± 8mim3( Jqp cos 20i - J^p sin 20i) , . 

aA<.(l + C) ' ^'"^ 
~± 4(m^ + m^) Jcp ±8m2m3(Jcpcos202 - Jcpsin202) 

23 ~ '^^^ ' ^ ' 

where a = S12C12S23C23. Clearly these quantities depend on Jcp, J^p in addition to the 
masses and Majorana phases. The coefficients for the RG evolution of Jcp and J^p are 
presented in Table llllt where the quantities Q^^c,s given by 

^ ^ ml + ml mj + mj 

Aml,^ Am2,^(l + C)' 
^ _ mims sin 201 m2m3 sin 202 /„o\ 
- AmlUl + Am^ ' ^^^^ 



atm 

^ _ mim3 cos 201 _^ m2m3 cos 2^^ 



Thus all the the quantities appearing in the evolution equations flTHl) - fl84l) have finite 
well-defined limits for ^13 — > in the Vj basis. 

Even if one starts with diagonal Ye (i.e. Ye = diag(|/e, y^, y^)) at the high scale, non-zero 

(n) („) 

off-diagonal elements of Ye will be generated through Eqs. ( |70l) - ( |73l) since Y^^Yj: is not 
diagonal. These off-diagonal elements will give additional contributions to the running of 
masses and mixing above and between the thresholds through F and «£• Since ag is fiavor 
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ReFsi 


^ zo 


-4c23Jcp/a 








-C23a/2 








ReF32 








1 














ImF2i 





-4s23</cp/a 





-S23«/2 





C23C12/S12 


-C23S12/C12 


ImFsi 





-4c23Jcp/a 





-C23a/2 





-S23C12/S12 


S23S12/C12 


ImF32 

















-1/(C23S23) 


-1/(C23S23) 


\BLE IV: 


Coefficients of Fjg 


in the RG evolution 


equations of all the ang 


;les (012, 0^3, 9: 



JcPjJcp and the Majorana phases (pi in the limit ^13 — > 0. The convention used here is a = 
•S12C12S23C23, and Jcp = (o/2)-si3cf3 sin 5. We neglect ye and compared to yr, and take vanishing 
flavor phases. 

diagonal, it will contribute to the running of ye, and yr, while off-diagonal conponents of 
F will contribute additional terms in the /^-functions of angles and phases, as tabulated in 
Table HVl These contributions will just get added to the Pfg contribution for the evolution of 
the quantities in Tables UT t llll t II VI Note that the Ffg coefficients are < C(l), whereas the Pfg 
coefficients are > 0{mf / Aml^^). Since the running is significant only when ^ Aml^^, 
in almost all the region of interest Pfg contributions dominate over the Ffg contribution. 

Note that the analytical expressions obtained in Eq. (1751) onwards, and those given in the 
tables, are valid only in the two extreme regions /i > M3 and /i < Mi. For the intermediate 

(n) 

energy scales, mj, will receive contributions from both k and Q. In the SM these two 
quantities have non-identical evolutions, as seen from Eqs. fl6T|) and fl62l) . and therefore the 
net evolution of Ye and is rather complicated. We perform it numerically in the next 
section. 
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V. ILLUSTRATIVE EXAMPLES OF RG RUNNING OF MASSES AND MIXING 



In this section we numerically calculate the RG evolution of the masses and mixing 
parameters within the Type-Ill seesaw model including the impact of running between the 
thresholds. This analysis is done by imposing suitable matching conditions fl25p at the 
thresholds. For illustration, we start at /iq = 10^^ GeV and choose the basis in which Ye is 
diagonal, so that f/pMNS = Uy- We further choose at this high scale to be the bimaximal 



mixing matrix t/,/,bimax [391, l69|], i.e. 612 = 623 = vr/4 and 6^13 = 0. This scenario is clearly 
inconsistent with the current data in the absence of RG evolution. We shall check if the 
radiative corrections to the masses and mixing angles can make it consistent with the data 
at the low scale. 

If the low energy theory in the complete energy range fi < fi^is the SM, then 612 decreases 
as the energy scale decreases, however the running is not sufficient to achieve compatibility 
with the low energy data. If the low energy theory is the MSSM, then 612 increases with 



decreasing energy sea! 



has been shown in 



59 



e jsG ]. so that compatibility with the data is not possible. However, it 



70 



71| in the context of Type-I seesaw mechanism, that the inclusion 
of threshold effects can make the mixing angle 612 decrease substantially as we go to lower 
energy scale and can give the correct values consistent with the Large Mixing Angle (LMA) 
solution. In this section we study the evolution from bi-maximal mixing at high scale in the 
context of Type-Ill seesaw scenario, including the seesaw threshold effects. 
We write the neutrino mass matrix as 

= ?7*bimaxdiag(mi, ms, m3)f/J t,.^„^ , (85) 

with 6e = Sfj, = 6r = a:t the high scale. Given the masses of the three fermion triplets and 
the light neutrino masses at the high scale, one can determine a at the high scale'^ that 
satisfies the seesaw relation = —(v'^/2)Y^'M^^Yy:. We then evolve the parameters using 
the analysis of Sec HVl 

Among the neutrino mixing angles, 612 is expected to be the most sensitive to RG effects. 
Table HI] shows that 612 is proportional to Qf2 and 1S12, which are in turn proportional to 



^ The solution for Is need not be unique, however any one of the solutions would suffice for the illustration. 
For practicality, we first choose an "trial" Ie, calculate the corresponding Ms from the seesaw relation, 
and then apply the basis transformation that makes Ms diagonal and takes the "trial" Is to its final 
form. 
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FIG. 3: RG evolution of mixing angles and mass squared differences, starting from bimaximal 
mixing at fiQ = 10^^ GeV, for normal mass ordering and hierarchical neutrino masses. The left 
panels represent the scenario where the Majorana phases vanish at /io- The right panel shows 
a representative case of nonzero Majorana phases {(pi = 89.0°, 02 = 0.4°) at ^q. The values of 
parameters at the high scale have been chosen such that the Am^'s and g2 at the low scale are 
reproduced. 

(m^/Amg^^i) as can be seen from Eqs. fl76l) and fl77p . For the other angles 6ij, the correspond- 
ing quantities Qfj and Sij are proportional to {mf/Aml^^), so the evolution of these angles 
is smaller. The direction of 612 evolution depends on the details of the Yukawa coupling 
matrix and masses of the heavy fermions. 

Since the values of Majorana phases at the low scale are completely unknown, we first 
consider the case where 01 = 02 = 0. In this case the CP violation will remain zero at all 
energy scales. The left panels of Fig. [3] show the running of mixing angles and mass squared 
differences for the normal mass ordering in this scenario. It is observed that ^12 in the 
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intermediate energy region changes more rapidly than in the extreme regions, however this 
change is in the opposite direction to what is required. As a result, bimaximal mixing at the 
high scale is not compatible with the low energy data in our model when the Majorana phases 
vanish. With nonzero Majorana phases, however, it is possible to achieve compatibility with 
the low scale data, as can be seen from the right panels of the figure. 

The lower panels of Fig. [3] show the evolution of mo, the lowest mass scale, and the two 
mass squared differences. As can be observed, the running of masses is quite substantial in 
Type-III seesaw, as compared to the SM, the MSSM js^, or the Type-I seesaw 59 1. Most of 
this running occurs in the intermediate energy range Mi < /i < M3, where threshold effects 
play a crucial role in enhancing the running. Note that the values of mo required to cause 
substantial running of mixing angles is quite small: in the case of vanishing (non vanishing) 
Majorana phases, we have taken mo = 0.04(0.01) eV at /i = fiQ. Thus, even at extremely 
small mo, substantial running of neutrino parameters can be present in the Type-III seesaw. 

The example of the bimaximal mixing discussed above was just for illustration. However, 
it brings out certain salient features of the RG running in Type-III seesaw scenario. The 
running of neutrino masses can be quite substantial here in the intermediate energy range. 
Moreover, threshold effects can enhance the extent of running of mixing angles, as well as 
the direction of the evolution, similar to the Type-I seesaw scenario [59|]. Majorana phases 
are also seen to play an important role in determining the extent and the direction of RG 
running of neutrino mixing parameters. 

In Fig.m we illustrate the RG evolution of parameters when the neutrino masses are quasi- 
degenerate. We have taken the parameter values at the high scale to achieve compatibility 
with the low scale data, without imposing any special symmetry. However in order to bring 
out certain salient features of the RG evolution that are independent of the threshold effects, 
we have chosen a small 613 value, |0i — ^2! ~ 7r/2, and Y^^Y-s to be almost diagonal in the 
charged lepton basis, with hierarchical eigenvalues. These conditions ensure that P21 and 
P31 are small, and S12 vanishes, so that from Table UTTl the evolution of (0i — 02) is extremely 
small. Thus |0i — 02| is expected to stay close to 7r/2 even after evolution, which is verified 
by the figure. Moreover, combined with mi ^ m2, the choice |0i — 02] ~ tt/2 makes Q12 
extremely small, thus restricting the ^12 evolution. 

It is observed that the running of 623 is now large, owing to ttiq/ Aml^^ ~ 1. This makes 
it possible to mimic maximal mixing accidentally, even if the mixing generated at the high 
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FIG. 4: RG evolution of mixing angles, mass squared differences, and CP violating phases, for 
quasi-degenerate neutrino masses and normal mass ordering. The values of parameters at the high 
scale have been chosen such that the Am^'s and 52 at the low scale are reproduced. Note that for 
the Majorana phases (pi, the regions (0° — 180°) and (180° — 360°) should be identified with each 
other. 

scale is arbitrary. The value of ^13 also quadruples from its high scale value. The Dirac 
phase, which was chosen to vanish at yUo, is generated by the RG evolution. The running of 
Dirac as well as Majorana phases is substantial between the thresholds. 

The right hand bottom panel of Fig. H] shows the evolution of rufs = ^/y]. \Uei\'^mf, 
the effective neutrino mass measured in the Tritium beta decay experiments 72], as well 
as rriee = \ Yli ^ei'^ilj effective neutrino Majorana mass in the neutrinoless double beta 
decay. Note that since 613 is small, mi ^ m2 ~ mg, and since |0i — 02| ~ 7r/2 in addition, we 
have rriee ~ mo cos 2^12. Also in the quasi-degenerate case, the sum of neutrino masses that 
is resticted by cosmology is ^ Sitlq. The large running of these masses suggests that, 

even if the beta decay experiments were to bound to < 0.3 eV, or the neutrinoless double 
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beta decay experiments were to bound niee to < 0.1 eV, or the cosmological observations 
were to restrict mo at the low scale to < 0.3 eV, the value of mo generated at the high scale 
can still be substantially larger. 

It is thus observed that in Type-III seesaw, the RG evolution of masses, angles as well as 
CP violating phases can be significant between the thresholds even at low mo values. The 
reason behind this, as well as the exact dependence of the evolution on the mass thresholds 
and Majorana phases, needs to be studied in further detail for a better understanding of the 
allowed neutrino parameter space at high energies. 

VI. SUMMARY AND CONCLUSIONS 

In this paper wc have studied the RG evolution of neutrino masses and mixing angles in 
the context of Type-Ill seesaw mechanism mediated by heavy fermions S transforming as 
triplets under SU(2) l- Tree level exchange of such particles gives rise to an effective operator 
i^blLlL<t>(f> below their lowest mass threshold. If one or more such triplets are present in the 
model, they affect the RG evolution of wavef unctions, masses and couplings. We compute 
these extra contributions using dimensional regularization and minimal subtraction scheme. 
We calculate the beta functions for the Yukawa couplings Ygj ^u-, and Y^i the SU(2)/, 
gauge coupling g2, the Higgs self-coupling A, the heavy fermion triplet mass matrix Ms, 
and finally the light neutrino mass matrix m,^. We do our calculation in the i?^ gauge and 
show the gauge invariance explicitly by demonstrating that the terms containing ^ are not 
present in the /3-functions. 

It is found that the presence of the triplets does not give rise to any additional diagram for 
the effective vertex k. However, the presence of these fields is felt indirectly in the running 
of K through their contribution to the evolution of the other quantities. Since the fermion 
triplets couple to W bosons, the evolution of the SU(2)i gauge coupling g2 is significantly 
affected, with more than two E triplets changing the sign of the (5 function for g2. This may 
also have implications for the unification of gauge couplings. In turn, the masses of the E's 
are also affected substantially due to the coupling with g2. 

We give the analytic expressions for the RG evolutions of the neutrino masses and 
mixing above the highest mass threshold and below the lowest one. We use a ba- 
sis Vj — {rrii, 9i2, 023,01^, (pi, Jcp, J'cp} instead of the commonly used basis Vs — 
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{mi,9i2,923,0i3,4>i,S}. The advantage of the Vj basis is that all the evolution equations 



5l|. 



are explicitly non-singular at all points in the parameter space including at ^13 = 

We consider the scenario with three triplets having non-degenerate masses and include 
the effect of successive decoupling of the heavy triplets at their respective mass thresholds by 
imposing suitable matching conditions at each threshold. We present illustrative examples 
of running of masses and mixings by numerical diagonalization of the effective neutrino mass 
matrix. Although the running of neutrino parameters is not very large in the SM, in our 
model the running can be large due to threshold effects of the heavy triplets. In particular we 
find that starting from bi-maximal mixing at a high scale it is possible to generate low scale 
values of masses and mixing angles for the normal hierarchical neutrino spectrum. However, 
this requires non-zero values of the Majorana phases. Indeed it is observed that threshold 
effects and Majorana phases can influence the evolution of the mixing angles significantly. 

We show that even in the case of hierarchical neutrinos, the RG evolution of neutrino 
masses and mixing between the thresholds can be substantial in the Type-Ill seesaw scenario. 
Moreover for quasi-degenerate neutrinos, the large running of masses implies that the value 
of mo at the high scale can be quite large, even if the mass related measurements from the 
beta decay, neutrinoless double beta decay, or cosmology, restrict its value at the low scale. 

In conclusion, this work studies threshold effects in the context of the Type-Ill seesaw 
mechanism. It is crucial for testing the viability of a high scale theory with low scale data. 
Indeed it is seen that theories that are excluded by the data in the absence of RG running 
can become viable once these effects are included. In order to determine the allowed neutrino 
parameter space at the high scale, a detailed exploration of the dependence of RG effects 
on various parameters is necessary. This is all the more important in view of the onset of 
the precision era in neutrino physics. 
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APPENDIX A: FEYNMAN RULES INVOLVING THE FERMION TRIPLET S 

In this appendix, we list the Feynmen rules involving the fermion triplets S. Following 
[3], we introduce the fermion flow arrow for the leptons, which is the gray arrow in the 
diagrams. The black arrows indicate the lepton number flow. However interactions involving 
S may violet lepton numbers and thus the S line does not carry any lepton flow arrow. 
For the lepton number conserving interactions, the two arrows are parallel for particles, 
and antiparallel for the charge-conjugate fields. The Feynman rules are also given for the 
effective operator in the low energy limit of the theory obtained by integrating out these 
heavy fermion triplets. 

1. Propagator 




2. 



Yukawa interactions 





S3' 
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Gauge boson interactions 
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5. Effective vertex k 



6. Counterterms for k 
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APPENDIX B: FEYNMAN RULES FOR THE SM FIELDS 



In this appendix, we list the Feynmen rules involving the SM fields only, also given in 
661], which are needed for our calculations. The directions of the arrows should be interpreted 
in the same way as stated at the beginning of Appendix |X1 
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1. Propagators 
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where C = G for B boson and ^ = ^2 for W boson. 



2. Yukawa interactions 





Similar Feynman rules, as those in the left panel, are there for Yukawa interactions of Ql-ur 
and QL-dn with the Higgs having coefficients Yu and Yd respectively. 
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3. Gauge boson — lepton interactions 




4. Gauge boson — Higgs interactions 



(f>a <Pa 
\ \ 

^ /W B X A* _ 

/ / 

06 0fe 



The vertices involving two Higgses and two gauge bosons are not shown since they do not 
appear exphcitly in our analysis. 



5. Higgs self-interaction 



0c 
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6. Counterterms 
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APPENDIX C: CALCULATION OF RENORMALIZATION CONSTANTS 

Here we show the Feynman diagrams contributing to the renormahzation constants of 
different quantities. Note that for particles in the loop, we suppress the flavor as well as the 
SU(2)i indices. 



1. Doublet Higgs wavefunction and mass (Z^ and (5m^) 
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Wavefunction and mass of fermion triplet {Zs and ^Me) 
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Right-handed charged lepton wavefunction (^efl) 
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6. lL'^(j) Yukawa vertex [Zyy) 
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7. The extra diagram contributing to Zx 
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Calculation of 
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